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1.  Brief  overview  of  technical  results 


This  grant  was  activated  on  March  1,  1990  with  the  project  period  of  three  years  ending 
on  February  28,  1993  (with  no-cost  extension  March  1,  1993  -  February  28,  1994).  The 
research  of  this  principal  investigator  has  continuously  been  supported  by  AFOSR  from  1980. 
During  that  entire  period  under  AFOSR  support,  the  principal  investigator  authored  or  coau¬ 
thored  212  publications,  among  them  5  book  contributions,  68  regular  journal  papers,  and  22 
conference  proceedings  papers;  the  rest  are  conference  papers. 

In  particular,  under  this  AFOSR  grant  #  90-0180,  more  than  65  new  papers  have  been 
published  by  this  principal  investigator  and  his  research  group,  among  them  15  papers  in  regu¬ 
lar  journals  [1-15],  4  book  contributions  [16-19],  8  conference  proceedings  papers  [20-27],  and 
32  conference  papers  [28-65]. 

Most  of  the  proposed  effects  are  novel  and  have  initiated  new  opportunities  in  the  field. 
The  work  by  this  principal  investigator  is  highly  credited  by  the  research  community  in  the  field. 
Within  the  1985-92,  for  example,  his  work  was  cited  for  more  than  620  times  (according  to  "Sci¬ 
ence  Citation  Index  ")  by  other  researchers  (only  those  papers  in  which  he  was  either  the  only 
or  the  first  author,  were  counted).  He  has  been  a  member  of  program  committees  and  a  panel 
member  of  several  technical  conferences  on  nonlinear  optics  and  quantum  electronics.  He  is 
an  editorial  board  member  of  two  international  journals  ("Optical  &  Acoustical  Review"  and 
"Nonlinear  Optical  Physics").  He  was  also  a  member  of  the  Natl.  Panel  on  "Math,  for  Material 
Sciences",  1991-92,  and  a  member  of  ASAT  US-USSR  Laser  Weapon  Verification  Meeting, 
June’91. 

Since  the  research  on  this  AFOSR  grant  started  on  March  1’90,  a  number  of  new  results 
were  obtained  by  this  principal  investigator  and  his  group  in  the  field  of  nonlinear  optics  and 
quantum  electronics.  The  research  progressed  basically  in  these  directions: 

l.i.  Pilot  theoretical  research  on  X-ray  nonlinear  optics,  including  resonant  saturation- 
related  effects  in  plasmas  (nonlinear  absorption  and  nonlinear  refractive  index),  and 
proposals  for  X-ray  third  harmonics  generation.  X-ray  laser  with  noncoherent  pumping, 
and  X-ray  four-wave  mixing. 

1  .ii.  Experimental  discovery  and  experimental  and  theoretical  research  on  dark  spatial  soli- 
tons. 

I.iii.  Theoretical  investigation  of  previously  discovered  by  this  PI  and  his  group  dispersion- 
related  multimode  amplification,  instabilities,  oscillations  and  chaos  in  nonlinear  coun- 
terpropagating  waves. 

I.iv.  Theoretical  research  on  the  fundamental  problem  of  gradient-field-induced  second- 
order  nonlinear  optical  processing  vacuum  due  to  the  photon-photon  scattering  of 
intense  laser  radiation  in  a  dc  magnetic  field. 


1  .V.  Collaborative  (with  Dr.  Khurgin’s  group)  theoretical  and  experimental  research  on  non¬ 
linear  optical  effects  (nonlinear  refractive  index,  anomalous  blue  shift  and  bistability)  in 
semiconductors. 

1.vi.  Other  research 


2.  Technical  reports  on  specific  projects 
2.i.  Research  on  X-Ray  Nonlinear  Optics 

In  our  research  under  this  AFOSR  grant,  we  explored  opportunities  for  X-ray  nonlinear 
optics  based  on  the  available  X-ray  lasers  as  potential  quasicoherent  sources  and  various  plas¬ 
mas  as  nonlinear  media  [7,  9,  10,  12,  14,  19,  23,  25,  26,  36,  46,  49,  50,  52-54,  56-58,  60-65]. 
The  proposed  effects  include  nonlinear  refractive  index  effects,  third-harmonic  generation, 
four-wave  mixing,  and  a  proposal  for  a  new  X-ray  laser  pumped  by  a  noncoherent  X-ray  pulse. 
This  work  was  done  by  P.  L.  Shkolnikov  and  this  P.  I.  The  work  on  X-ray  four-wave  mixing  was 
done  in  collaboration  with  P.  L  Hagelstein  and  M.  H.  Muendel  of  MIT. 

After  the  invention  of  optical  lasers,  nonlinear  optics  of  visible  domain  has  developed  into 
one  of  the  most  important  fields  of  optics  with  numerous  applications  to  both  science  and  tech¬ 
nology.  Now  its  shorter-wavelength  counterpart,  X-ray  nonlinear  optics  with  similarly  great 
potentials  and  opportunities,  is  making  its  first  steps.  Although  the  very  first  publications  on  X- 
ray  nonlinear  optics  (primarily  about  the  parametric  conversion  of  X-rays  in  solids)  appeared  in 
late  60’s-early  70’s,  the  entire  field  of  X-ray  nonlinear  optics  had  to  wait  for  X-ray  lasers  (XRLs) 
to  come. 

The  first  X-ray  amplification  was  observed  in  1984  at  LLNL  and  PPPL.  By  now,  more 
than  fifty  X-ray  laser  lines  have  been  reported,  with  the  wavelengths  ranging  from  326.5  % 
down  to  35.6  ?l,  with  the  output  power  up  to  several  MW.  Such  rapid  progress  in  XRL  research 
revived  interest  in  X-ray  nonlinear  optics.  Apart  from  pure  scientific  interest  (in  particular,  for 
X-ray  laser  spectroscopy).  X-ray  nonlinear  optics  also  appears  to  have  a  great  potential  for 
applications  to  X-ray  laser  technology.  A  few  X-ray  nonlinear  effects  were  preliminary  con¬ 
sidered  in  1989  and  1990.  A  number  of  researchers  working  on  X-ray  lasers  mentioned  X-ray 
nonlinear  optics  as  one  of  the  most  important  future  applications  of  X-ray  lasers. 

This  principal  investigator  was  among  the  first  to  realize  importance  and  timeliness  of  X- 
ray  nonlinear  optics,  and  the  first  to  propose  full-scale  research  on  the  subject.  Under  this 
AFOSR  grant,  19  papers  have  been  published  on  the  subject,  among  them  5  papers  in  regular 
journals  [7,  9,  10,  12,  14],  a  book  contribution  [19],  two  conference  proceedings  paper  [23,  25], 
and  11  conference  papers  (among  them  two  invited  [57,58]).  Together  with  the  early  work, 
these  results  mark  the  birth  of  the  theoretical  X-ray  nonlinear  optics. 


From  the  very  beginning  we  realized  that  the  main  issue  for  X-ray  nonlinear  optics  is  not 
the  theory  of  nonlinear  propagation  per  se  (which  in  its  major  features  is  substantially  similar  to 
that  of  visible  domain),  but  the  search  for  media  in  which  X-ray  nonlinear  effects  may  be 
observed.  Our  very  first  estimations  of  X-ray  nonlinearities  showed  that  the  main  hope  to 
achieve  feasible  effects  at  such  short  wavelengths,  is  to  use  resonances.  We  believe  now  that 
the  entire  business  of  X-ray  nonlinear  optics  is  about  resonant  couples,  i.  e.  couples  of  XRL  line 
-f-  resonant  transitions  in  a  nonlinear  media.  From  this  point  of  view,  plasmas  turned  out  to  be 
the  first  candidates  to  be  considered  as  X-ray  nonlinear  media.  That  is  because  there  are 
much  larger  selection  of  ionic  spectra  than  that  of  neutral  atoms.  Indeed,  Z-1  different  degrees 
of  ionization  (that  is,  Z-1  different  ionic  spectra)  are  possible  for  an  atom  with  the  atomic 
number  (nuclear  charge)  Z.  Thus,  if  the  total  number  of  elements  is  Zmax.  the  total  number  of 
ionic  spectra  is  ~Z^ax^.  '•  e-  'Tiore  than  5000.  Moreover,  majority  of  ionic  spectra  consists 
predominantly  of  X-ray  lines  related  to  what  may  be  called  “Z-scaling".  For  example,  the  spec¬ 
trum  of  a  hydrogen-like  ion  (i.  e.  an  ion  with  only  one  electron)  with  nuclear  charge  Z  is  roughly 
similar  to  the  spectrum  of  the  atomic  hydrogen.  The  major  difference  is  that  all  energies  are 
scaled  up  by  the  factor  of  Z^ . 

Thus,  until  now,  our  research  interests  have  been  concentrated  on  identification  and  esti¬ 
mation  of  couples  "X-ray  laser  -i-  plasma"  resonant  enough  for  the  feasibility  of  various  X-ray 
nonlinear  effects.  In  the  framework  of  our  theoretical  research  on  X-ray  nonlinear  optics,  we 
have  identified  resonant  nonlinear  media  (mainly  plasmas)  for  about  1/3  of  existing  XRL  lines. 
For  these  resonant  media,  we  have  also  theoretically  demonstrated  the  feasibility  of  several 
basic  X-ray  nonlinear  optical  effects  such  as 

(i)  nonlinear  refraction  and  absorption  saturation, 

(ii)  third-harmonic  generation,  and 

(iii)  four-wave  mixing. 

We  believe  that  our  results  are  promising  enough  to  encourage  experiments  on  X-ray 
nonlinear  optics.  We  have  also  proposed 

(iv)  a  new  X-ray  laser  in  Li  vapor  pumped  by  a  short  noncoherent  X-ray  pulse. 

2.i.l.  Saturation-related  X-ray  nonlinear  effects  in  plasma. 

Saturation-related  X-ray  resonant  nonlinear  effects  in  plasma  were  chosen  by  us  to  be 
considered  first.  The  reason  for  that  was  that  we  wanted  to  use  X-ray  resonant  nonlinearities 
inside  an  XRL  itself.  Our  results  suggested  [7]  that  significant  nonlinear  change  is  feasible. 
For  example,  the  output  intensity  of  Se  XRL  line  of  206  %  is  much  larger  that  the  saturation 
intensity  for  the  laser  active  medium.  It  may  lead  to  an  order  of  magnitude  decrease  in  gain 
coefficient  due  to  saturation.  Furthermore,  "nonlinear  distance"  turns  out  to  be  comparable  to 
the  length  of  the  XRL  active  medium.  [Here  the  "nonlinear  distance"  is  the  distance  at  which 
significant  change  in  the  phase  front  (up  to  2  ti  across  the  beam)  may  result  from  nonlinear 
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refraction]. 

Transition  configuration  used  in  XRL’s  is  not  very  suitable  when  one  wants  to  seek  exam¬ 
ples  of  coherent  X-rays  with  separately  prepared  plasmas.  Indeed,  in  the  case  of  laser  active 
medium,  both  of  resonant  levels  are  excited  levels,  and  the  lower  level  has  to  be  populated,  as 
it  happens  inside  XRL’s.  Outside  of  an  XRL’s,  both  of  them  would  be  empty.  Thus,  to  study 
nonlinear  optics  interactions  of  XRL  radiation,  it  would  be  much  more  desirable  to  find  much 
cooler  plasmas  ,  and  with  resonant  ionic  transitions  from  their  ground  (not  excited)  levels. 

A  plasma  is  a  good  candidate  for  saturation-related  X-ray  nonlinear  effects  if  the  resonant 
detuning  |v-Vol  between  XRL  frequency  v  and  frequency  Vq  of  some  transition  from  the 
ground  to  an  excited  level  of  the  plasma  ions,  is  smaller  than  the  Doppler  full  width  on  half 
maximum  of  the  respective  plasma  transition  Av^.  Search  for  such  resonances  is  complicated 
by  the  lack  of  atomic  information  (wavelengths,  oscillator  and  collisional  strengths),  especially 
on  highly-ionized  atoms.  We  have  identified  such  resonant  couples  "  XRL  line  +  ion  "  for  about 
one  third  of  all  the  reported  XRL  lines  (Table  1).  As  an  example,  we  present  our  results  for 
two  X-ray  lasers:  Se  and  Ge  (Table  2).  The  last  column  in  the  Table  2  contains  the  estimates  of 
the  "nonlinear  length"  for  given  intensity  /  =  5x/^  where  4  is  the  saturation  inten¬ 

sity.  One  can  readily  see  from  Table  2  that  nonlinear  effects  resulting  from  nonlinear  refractive 
index  can  be  observed  at  the  Ge  and  Se  XRL  intensities  even  much  lower  than  the  intensities 
available  now.  The  expected  effects  based  on  nonlinear  refractive  index  are:  self-focusing  and 
"running  foci"  (although,  because  of  significant  probability  of  multiphoton  processes,  we  do  not 
expect  a  pronounced  wave  collapse,  in  contrast  with  the  situation  in  liquids  and  solids):  self¬ 
trapping  of  XRL  radiation,  self-bending  and  four-wave  mixing.  These  effects  may  also  be  instru¬ 
mental  in  direct  measurement  of  nonlinear  refractive  index,  plasma  diagnostics,  and  for  phase- 
conjugation  amplification.  It  may  also  be  possible  that  self-  transparency  and  related  2k- 
solitons  can  be  attained  similar  to  that  in  the  optical  domain.  The  experimental  observation  of 
XRNE’s  in  resonant  couples  of  Table  2  is  greatly  facilitated  by  the  fact  that  the  required  param¬ 
eters  Te,Ne,  and  L  of  plasmas  can  be  attained  by  using  standard  discharge  devices. 
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XRL  Xm 

Plasma  XqI^ 

A{s~'') 

Eton  {eV) 

1. 

Ge^'^  236.26 

Ar  XIII  236.27 

3x10® 

686 

2. 

5^24+  220.28 

Sc  IV  220.280 

6x10® 

73 

3. 

5^24+  209.78 

Fe  X  209.776 

1.3x10® 

262 

4. 

209.78 

Cl  XIII  209.81 

10^° 

656 

5. 

196.06 

Fe  VII  196.046 

X 

o 

125 

6. 

Ge^'^  196.06 

Na  III  196.054 

bo 

X 

O 

CO 

72 

7. 

182.173 

Na  IV  182.123 

1.4x  10^® 

99 

8. 

105.69 

Na  IV  105.6867 

2.4x10^° 

99 

9. 

105.69 

CrVIII  105.69 

4.3x10^® 

185 

10. 

0^+  102.355 

Ni  IX  102.340 

CO 

X 

193 

11. 

71.00 

Mg  VIII  71 .007 

- 

266 

12. 

44.83 

Si  X  44.83 

- 

401 

13. 

W4&+  43.18 

C/  7X43.168 

5x  10^^ 

400 

Table  1 .  Resonant  couples  "XRL  radiation  -  plasma-ion  transition".  A  denotes  the  correspond¬ 
ing  transition  probability,  Eio„  -  the  plasma  ionization  potential,  X  and  Xq  -  the 
wavelengths. 


Couple 

/ 

(W/cm^) 

r 

1. 

10® 

600 

3x10“^ 

6 

4. 

2x10® 

5x10® 

o 

1 

18 

Table  2.  Saturation-related  X-ray  nonlinear  effects  in  plasmas. 

2.i.2.  X-ray  third-harmonic  generation  in  plasmas  of  alkali-like  ions 

Absorption  saturation  and  nonlinear  refraction  were  our  first  research  objects  in  X-ray 
nonlinear  optics.  Another,  probably  even  more  important  effect,  is  higher-harmonic  generation, 
since  it  is  a  method  of  choice  to  generate  shorter-wavelength  coherent  radiation.  Therefore, 
our  next  problem  was  to  identify  resonant  media  for  X-ray  third-harmonic  generation  [10]. 
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(Since  plasma  is  isotropic,  the  third  harmonic  is  the  lowest  possible  higher-harmonic,  in  the 
first-order  approximation  in  velocities  of  free  electrons  and  for  not  very  strong  optical  field). 

We  face  here  the  same  problem  of  resonant  coupling.  But  in  this  case  the  complication  is 
in  a  sense  tripled.  Indeed,  it  is  well  known  that  the  best  resonant  conditions  for  THG  are  the 
closest  possible  two-photon  resonance  and  near  (to  a  few  linewidths,  to  avoid  strong  absorp¬ 
tion)  one-  and  three-photon  resonances  Strong  optical  coupling  between  levels  of  interest  is 
also  desirable.  Because  there  is  no  tunable  XRL  available,  one  has  to  seek  for  plasmas  that 
may  allow  for  resonances  to  existing  XRL  lines.  This  search  is  complicated  by  insufficient 
information  on  energy  levels  and  transition  probabilities  of  ions.  In  this  respect,  plasmas  of 
alkali-like  ions  present  perhaps  the  best  opportunity  because  of  readily  available  atomic  data. 
We  have  identified  seven  resonant  couples  for  X-ray  third-harmonic  generation  with  plasmas 
on  alkali-like  ions  (see  Table  3). 

To  estimate  conversion  efficiency,  we  use  standard  results  of  the  theory  of  THG  in  gasses 
and  vapors,  modifying  them  mainly  to  account  for  the  specifics  of  the  plasma  dispersion.  The 
X-ray  refractive  index  due  to  free  plasma  electrons  is  np(3v)  =  1  -to^/co^  where  co  is  the  (angu¬ 
lar)  frequency  of  the  incident  XRL  radiation,  cOp  =  (iV^e^/eom)^^  is  the  plasma  frequency,  and 
Ng  is  the  plasma  electron  density.  For  all  but  one  the  couples  "XRL  line-plasma"  in  Table  3,  this 
free-electron  component  of  the  refractive  index  is  much  larger  than  the  resonant  (bound- 
electron)  one,  which  can  therefore  be  neglected.  In  our  estimates  we  also  assumed  that  the 
plasma  is  homogeneous  enough  for  us  to  neglect  the  refraction  due  to  inhomogeneous  elec¬ 
tron  distribution  across  the  laser  beam.  As  a  result,  the  phase  mismatch 
AA:p  =  67t?."^[n/,(3(o)-«p(a))]  where  X  denotes  the  XRL  radiation  wavelength,  is  positive.  The 
results  of  our  calculations  can  be  found  in  Table  3  in  terms  of  the  XRL  intensities  required  to 
attain  Cgff^  1 (collimated  beams  assumed).  They  show  that  the  observable  THG  intensities 
can  be  attained  by  available  XRL’s. 
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Resonant  couples 
XRL  line  X{nm) 
plasma /Eioni^V) 

(X/3)(nm) 

for 

Buffer  ions 
a/E'UeV) 

I{Wcm-^) 

for 

28.646 

KIXAie 

9.5487 

8x10^® 

- 

- 

20.638 

ScXI/250 

6.8793 

2x10^^ 

- 

- 

6.0810 

5x10^2 

Si  VII 

0.06/247 

2x10^^* 

18.2097 

Ne  W///239 

6.0699 

2x10^® 

Be  IV 

0.20/218 

10^° 

5.167 

2x1 0^"^ 

- 

- 

3.31 

5x10^^ 

- 

- 

Em®®-"  6.583 

Ga  XXI  mi 

2.194 

2x10^^ 

TiXIV 

0.06/863 

7x10^° 

Na^^  5.4194 

Cl  XV /809 

1 .8065 

'io'® 

- 

- 

Table  3.  X-ray  THG  in  plasmas,  a  is  the  ratio  of  the  required  buffer  ion  density  to  the  electron 
density,  and  and  are  the  ionization  potentials  of  the  "main"  and  "buffer" 
plasmas,  respectively.  The  asterisk  marks  the  result  obtained  for  ion  density 
10^®  cm -2. 


(For  example,  the  power  of  XRL  required  for  10~®  in  the  loose  focusing  limit  can 

be  estimated  as  several  MW  which  is  close  to  the  peak  power  of  the  existing  XRL.  It 
corresponds  to  the  input  energy  of  several  hundreds  of  microjoules  in  a  100  ps  pulse). 

One  of  the  ways  to  improve  the  X-ray  THG  is  to  facilitate  ideal  phase  matching  condition 
M  =  0.  Similarly  to  the  optical  THG,  this  can  be  done  by  adding  a  "buffer"  medium  or,  in  our 
case,  buffer  ions.  Let  some  transition  from  the  ground  level  of  the  buffer  ions  be  resonant  to 
the  third  harmonic,  X,/3  <  Xo  where  Xq  is  the  central  wavelength  of  the  transition,  and  let  no 
transition  be  in  close  resonance  to  the  fundamental  harmonic.  Then  the  resonant  refraction  by 
buffer  ions  may  compensate  for  free-electron  phase  mismatch  if 

1  -nto^3co)  =  2x10-®(X/3)(//Av^)N;  w{x) 

with  the  buffer  ion  density  N'i  in  cm"®  and  X.  in  cm.  Required  plasma  conditions  may  be 


attained  in  some  discharge  devices.  This  way  one  can  achieve  several  orders  of  magnitude 
enchancement.  (For  example,  for  C  XRL,  Ne  VIII  plasma  as  a  nonlinear  medium  and  Be  IV  as 
buffer  ions  the  intensity  required  to  attain  Qj9-=10"®  is  four  orders  of  magnitude  smaller  com¬ 
pared  to  unbuffered  case.) 

Our  quantitative  results  should  be  viewed  as  preliminary  estimates,  in  particular  because 
at  this  stage  some  of  the  competing  processes  have  been  neglected.  Nevertheless,  we  expect 
substantial  third-harmonic  intensities  at  the  state  of  the  art  of  both  X-ray  lasers  and  plasma 
technology.  Further  significant  enhancement  of  the  efficiency  (up  to  several  orders  of  magni¬ 
tude  for  the  same  XRL  intensity)  require  search  for  better  resonant  couples  "XRL  line-plasma". 

2.i.3.  X-ray  laser  frequency  near-doubling  and  generation 
of  tunable  coherent  X-rays  by  four- wave  mixing  in  plasma 

Another  third-order  nonlinear  effect,  which  may  lead  to  X-ray  laser  frequency  upconver- 
sion,  is  four-wave  mixing  of  coherent  X-ray  and  optical  radiation  [14].  One  could  expect  high 
conversion  efficiency  for  this  process  because  (i)  powerful  optical  lasers  may  be  used;  (ii)  one 
of  the  participating  frequencies  is  relatively  low  (optical):  (iii)  close  resonances  are  more  likely 
found.  Moreover,  this  process  may  appear  instrumental  in  generating  tunable  coherent  X-rays. 

We  identify  and  evaluate  plasmas,  in  which  efficient  X-ray  laser  frequency  near-doubling  is 
expected  for  a  number  of  available  X-ray  lasers  by  means  of  resonant  four-wave  mixing  with 
Nd  or  KrF  lasers.  In  some  of  these  plasma,  four-wave  mixing  of  coherent  X-rays  and  tunable 
optical  radiation  may  result  in  tunable  coherent  X-ray  radiation  powerful  sufficient  for  X-ray 
laser  spectroscopy. 

Of  the  various  four-wave  mixing  processes,  we  have  chosen  to  study  difference-frequency 
mixing  since  in  this  case  it  seems  relatively  easy  to  attain  optimal  phase  matching.  As  holds 
for  X-ray  nonlinear  optics  in  general,  only  resonantly-enchanced  effects  may  be  observed  at 
the  present  level  of  X-ray  laser  technology.  We  are  therefore  interested  in  resonant  nonlinear 
media  (plasmas)  for  the  following  difference-frequency  mixing  process; 

(0  =  2(i)xRL-(>iopt  (A) 

Here  (Hxrl  is  the  frequency  of  an  XRL,  and  cOop,  is  the  frequency  of  an  optical  laser.  For  the 
purpose  of  XRL  frequency  near-doubling,  we  consider  processes  (A)  with  two  particular  optical 
lasers  :  a  powerful  Nd  laser,  which  is  used  to  pump  most  XRL’s,  or  a  very  bright,  short-pulse 
KrF  laser.  (In  the  latter  case,  the  converted  X-ray  pulse  duration  may  be  much  shorter  than  the 
duration  of  the  incident  X-ray  pulse.)  If  a  tunable  optical  laser  is  used  instead,  the  result  is  tun¬ 
able  coherent  X-ray  radiation  at  almost  doubled  XRL  frequency. 

Ions  of  two  isoelectronic  sequences  could  be  most  useful  for  X-ray  four-wave  mixing.  For 
many  C-like  ions,  the  energies  of  2s^2p^  -  2s2p^  and  2s2p^  ^P°  -  ^P  transi¬ 

tions  are  very  close  to  each  other,  providing  one  with  one-  and  two-X-ray-photon  resonances. 


Transitions  from  2p'^  levels  to  the  nearest  lower  level  2s 2p^  ^S-\  correspond  to  X-ray  radi¬ 
ation  if  the  previous  two  do,  and  therefore  no  three-photon  resonances  are  possible  for 
processes  (A)  in  C-like  ions.  By  the  same  token,  tunable  X-ray  radiation  may  be  generated  in 
these  ions  through  process  (A)  for  a  broad  range  of  optical  frequencies  with  good  efficiency. 

Na-like  ions  may  provide  all  three  resonances  for  processes  (A).  Also,  the  plasma  tem¬ 
perature  required  for  substantial  presence  of  the  Na-like  ionization  stage  is  relatively  low  due  to 
low  ionization  potentials  of  Na-like  ions.  A  disadvantage  of  Na-like  ions  as  X-ray  four-wave 
mixing  nonlinear  media  is  that  initial  levels  of  nonlinear  transitions  should  be  excited  levels 
(only  then  are  simultaneous  one-  and  two-photon  resonances  possible),  and  it  may  be  difficult 
to  populate  them  significantly. 

Resonant  combinations  of  lasers  and  ions,  with  estimates  of  conversion  efficiency,  are 
listed  in  Table  4.  The  tight-focusing  limit  ^  <*:L  is  assumed  {b  is  the  confocal  parameter  and  L 
is  the  length  of  a  plasma  cell).  To  evaluate  the  minimum  plasma  length  required,  we  make  use 
of  the  feasible  spot  size  w  =  2  p/n  of  the  155?lXRL  (confocal  parameter  =  0.15  cm).  Simi¬ 
lar  values  of  the  confocal  parameter  are  readily  attainable  for  the  Nd  laser  fundamental  har¬ 
monic.  We  assume  b  =  0A5cm  for  all  the  participating  beams.  It  means  that  plasma  length 
Lrnin  =  1 0  Z?  =  1 .5  cm  satisfies  the  tight-focusing  condition. 


Medium  Eio^ 

Lasers 

K 

X 

^eff 

\.Ar  VIII  143  eV 

236.26%-Nd  1.06  pm 

3x10”“ 

11 9.46  a 

it 

2.KIX  176  eV 

Se^‘^  209.78%- Nd  0.53 |im 

10-52 

1 06.99 

•* 

3.  Ca  X  271  eV 

Se^^  206.38%- Nd  0.53pm 

5x10-5^ 

105.24^ 

* 

4.K  IX  176  cV 

5^2^  206.38%- Nd  0.266pm 

10-51 

1 07.35 

5.K  IX  176  eV 

182.173^- AW  0.263pm 

3x10-55 

94.29  a 

★ 

6.  Ca  X  271  eV 

1 55.0 0.5320  pm 

10-54 

78.65  /I 

it 

7.  y  XIII  336  eV 

99.36%- Nd  1.06  pm 

10-53 

49.91  ^ 

★ 

8.  Cu  XIX  670  eV 

44.83  AW  1.06  pm 

4x10-55 

22.46 

0.01 

9.  Ca  XV  894  eV 

182.097A-A:rF  0.2484  pm 

10-53 

94.51 

* 

lO.Ar  X7//  686  eV 

Se^"^  209.78%- KrF  0.2484pm 

6x10-57 

109.51 

0.03 

Table  4.  X-ray  laser  frequency  near-doubling  by  four-wave  mixing  with  optical  laser  in  plasma. 
*  denotes  conversion  efficiency  near  1  (with  competing  processes  neglected). 


In  order  to  improve  phase  matching,  the  product  \b  Ak  \  should  be  made  as  small  as  possible. 
For  collinear  beams  with  fixed  confocal  parameters,  this  can  be  done  by  adjusting  of  the  elec¬ 
tron  density.  For  example,  for  assumed  b=0A5cm  and  X2  =  2.484x1 0“^  cm  (KrF  laser), 

=  for  =  cm~^.  It  yields  |F2l^=4  which  is  only  two  times  smaller  that  the 

optimal  value  of  8  (for  zero  phase  mismatch  M  =  0). 

For  the  sake  of  uniformity,  all  the  values  of  Cgff  listed  in  Table  4  are  estimated  for 
Ne  =  ^0^^  cm~^  and  1^2!^  =4.  Our  quantitative  results  should  be  viewed  as  only  order-of- 
magnitude  estimates,  largely  because  at  this  stage  we  neglected  all  the  competing  processes. 
(It  is  worth  noting  ,  however,  that  substantial  photoabsorption  is  unlikely  to  occur  for  all  the 
cases  considered,  due  to  relatively  low  density  and  small  length  of  plasma,  and  due  to  the 
absence  of  close  resonances.)  Nonetheless,  these  estimates  allow  one  to  expect  substantial 
output  power  at  almost  doubled  frequency  for  a  number  of  available  XRL’s. 

Tunable  X-ray  radiation  may  result  from  four-wave  mixing  (A)  with  tunable  optical  lasers. 
In  particular,  a  tunable  laser  with  the  wavelength  near  5000  ^ and  the  power  P2  would  produce 
broadly  tunable  output  P  =  10"®xF2  at  almost  doubled  XRL  frequency  by  mixing  with  XRL 
in  Ca  XV  plasma  or  with  209.78 XRL  in  Ar  XIII  plasma  (from  lines  9  and  10  of  Table  4, 
respectively).  For  P2  of  a  few  MW,  converted  power  of  tens  of  watts  which  may  be  sufficient 
for  the  purposes  of  X-ray  linear  spectroscopy. 


2.i.4.  Proposal  for  19.9  nm  laser  in  Li  pumped  by  a  noncoherent  X-ray  pulse. 


All  X-ray  lasers  reported  to  date  operate  in  hot  and  dense  plasma  which  requires  very 
high  pumping  power  and  results  in  broad  (-lO”^)  bandwidth  of  amplified  X-rays.  A  variety  of 
pumping  schemes  have  been  proposed  for  the  last  three  decades  to  attain  X-ray  lasing  in  a 
cold  medium.  Most  of  them  rely  on  removal  of  atomic  inner  electrons  by  either  noncoherent 
X-rays  or  by  electron  impact.  An  X-ray  laser  pumped  through  inner-shell  photoionization,  how¬ 
ever,  has  not  been  realized  yet.  One  of  the  main  obstacles  to  proposed  lasing  to  the  ionic 
ground  level  [49]  has  been  that  such  pumping  would  create  free  electrons  energetic  and 
numerous  enough  to  destroy  population  inversion  by  electron-impact  ionization  of  neutral 
atoms. 

We  have  theoretically  demonstrated  [12]  how  the  problem  of  19.9  nm  lasing  on  the 
^s2p-^s^  transition  in  pumped  by  the  inner-shell  photoionization  and  photoexcitation  on 
neutral  Li  atoms,  can  be  solved.  To  limit  the  detrimental  effect  of  the  photoelectrons  created 
by  the  pumping  radiation,  on  the  population  inversion,  we  suggest  making  use  of  a  short  soft- 
X-ray  pulse  with  a  narrow-band  spectrum  centered  at  the  Li  [I5  {3s  3p)  resonance  to 

excite  ground-state  Li  atoms.  Subsequent  fast  autoionization  of  part  of  the  inner-shell  excited 
Li  atoms  to  the  upper  laser  level  {Li^ ‘\s^2p '' P)  is  complemented  by  the  direct  1s- 
photoionization  of  Li  atoms  to  the  upper  laser  level.  Suggested  pumping  scheme  relies  upon 
the  relatively  large  both  photoabsorptTon  cross-section  and  the  probability  of  the  inner-shell 
(auto)ionization  to  the  ^s2p  '' P  state  in  the  vicinity  of  g  resonance  of  neutral  Li  (  =71  eV 
above  the  ground  level).  At  the  same  time,  our  calculations  show  that  the  energy  distribution  of 
emerging  photoelectrons  does  not  allow  them  to  significantly  destroy  the  population  inversion. 
The  spectral  intensity  of  the  short  pumping  pulse  required  to  attain  significant  gain  seems 
feasible  at  the  the  state  of  the  art  of  the  laser-plasma  technology.  It  is  also  possible  that  much 
lower  X-ray  pulse  spectral  intensity  may  be  required  to  pump  the  excited  Li  atoms  by  a  broad¬ 
band  short  X-ray  pulse. 

It  is  well  known  that  the  absorption  of  soft  X-ray  radiation  (the  photon  energy  in  the  region 
57-100  eV)  by  Li  atoms  may  produce  double-excited  states  with  energies  much  higher  than  the 
2s-ionization  potential  of  5.39  eV  of  neutral  Li.  A  Li  atom  in  many  of  these  states  is  unstable 
against  autoionization:  within  a  very  short  time  (  10-100  ps  )  it  emits  a  free  electron  whereby 
turning  into  a  ion  in  one  of  the  excited  states. 

Double-excited  autoionizing  Li  atomic  states  correspond  to  broad  resonances  in  the  pho¬ 
toabsorption  spectra  of  Li  vapor  for  the  energy  of  the  incident  photons  above  58  eV.  For  the 
incident  photon  energy  large  enough  to  ionize  a  ground-state  Li  atom  to  the  Li^  ^s2p  ^P 
State,  the  strongest  such  a  resonance  corresponds  to  the  [l5'(3s3p)  P]  P  (71.14  eV 
above  the  Li  ground  level)  state  of  Li  atoms.  The  energy  interval  between  P' =71.05  and 
P"  =71 .25eV  approximately  corresponds  to  the  width  of  this  resonance.  Absorption  of  a  pho¬ 
ton  with  such  energy  would  immediately  result  in  either  (i)  excitation  of  Li  atom  to  the 
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[^s  {3s3p)^P]^P°  state  followed  by  the  autoionization  to  one  of  the  four  ionic  state: 
1^25  ^5,  1525  ^S,  ^s2p  ^P,  and  l5  2p^P  ,  and  by  freeing  of  an  electron  with  the  energy  of 
approximately  5.0  eV,  6.7  eV,  3.5  eV,  or  4.5  eV,  respectively:  or  (ii)  direct  (nonresonant)  inner- 
shell  ionization  to  the  same  ionic  and  photoelectron  states;  or  (iii)  2s-ionization  i.  e.  creation  of 
a  ion  in  the  ground  state  and  a  free  65.7  eV-electron. 

The  process  (iii)  creates  ground-state  Li  ions  whereby  decreasing  the  population  inver¬ 
sion.  However,  the  2s-photoionization  cross  section  is  very  low  for  the  photon  energy  so  high 
above  the  2s-ionization  potential  of  5.39  eV.  As  a  result,  ground  level  population  due  to 
the  direct  1  s-photoionization  does  not  exceed  0.013  of  the  population  of  the  ionic  excited  levels 
and  may  be  neglected. 

Processes  (i)-(iii)  are  directly  initiated  by  the  pumping  radiation.  In  turn,  they  bring  about  a 
variety  of  secondary  processes.  We  have  evaluated  the  most  important  of  them:  a)  the  2s- 
ionization  of  Li  atoms  by  electron  impact  (which  has  been  considered  as  the  main  danger  to 
the  population  inversion);  b)  collisional  and  optical  transitions  from  the  upper  laser  level  to  lev¬ 
els  other  than  the  lower  laser  level;  c)  the  electron-impact  excitation  of  the  ground-state  Li 
atoms.  Out  calculations  show  that  all  these  processes  would  hardly  diminish  the  inversion 
significantly:  before  substantial  amount  of  Li  atoms  could  be  ionized  by  photoelectrons,  almost 
all  of  these  electrons  would  lose  their  energy  through  the  collisional  excitation  of  the  ground- 
state  Li  atoms  to  the  l5^2p  level  (process  c)),  due  to  very  large  cross  section  of  this  process. 

Estimated  pumping  intensity  seems  moderate.  For  instance,  incident  noncoherent  X-ray 
spectral  intensity  required  to  attain  gain-length  product  of  =  1 0  in  a  traveling-wave  design,  is 
radiated  by  105  eV  black  body.  Such  a  temperature  is  readily  attainable  in  laser-produced 
plasma  (one  may  compare  this  temperature  with  - 1500  eV  temperature  of  existing  X-ray  laser 
active  media). 

For  considered  pumping  by  the  1  s-photoionization  of  Li  atoms  from  their  ground  state 
15^25,  the  probability  of  the  ion  to  be  left  in  the  1s2p  configuration  is  significantly  smaller  than 
in  the  1s2s  configuration.  Only  in  the  vicinity  of  [l^  {35 3p)  ^P  ]  ^P°  resonance,  this  proportion 
changes  somewhat  to  favor  19.9  nm  lasing.  As  a  result,  such  lasing  would  require  narrow- 
band  pumping.  However,  some  recent  experimental  results  on  the  soft-X-ray  inner-shell  ioni¬ 
zation  of  another  alkali,  Na,  allow  us  to  suggest  that  the  proposed  19.9  nm  X-ray  laser  in  Li 
may  be  efficiently  pumped  by  photoionization  of  excited  l5^2p  Li  atoms  with  a  short  X-ray 
pulse  of  a  much  broader  spectrum.  We  are  going  to  consider  such  pumping  in  more  detail 
later. 


2.i.5.  New  subjects  in  X-ray  nonlinear  optics 

Close  to  the  end  of  this  grant  period,  we  began  research  along  a  new  line  in  X-ray  non¬ 
linear  optics,  with  the  intent  to  explore  it  in  detail  under  our  continued  (current  at  this  point) 
AFOSR  grant.  This  new  research  is  related  to  the  recently  discovered  nonlinear  effect  -  very 
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high  order  harmonic  generation  of  optical  lasers,  with  the  highest  harmonics  reaching  into  soft- 
X-ray  domain.  We  started  exploring  new  nonperturbative  mechanism  of  this  process,  the  so 
called  super-driven  two-level  atom  [27,60,64,65].  We  also  began  research  on  optimal  phase¬ 
matching  conditions  for  frequency  upconversion  of  optical  laser  radiation  into  X-ray  domain  in 
plasma  [63].  Our  detailed  results  on  this  new  subjects  which  were  obtained  under  current  (at 
this  point)  AFOSR  grant  have  been  published  in  a  few  journal  papers. 

2.ii.  Spatial  Dark  Solitons  in  Self-Defocusing  Materials 

In  our  research,  a  new  phenomenon  of  spatial  dark  solitons  (SDS)  have  been  experimentally 
discovered  and  theoretically  explained  [1,4,24,  31,  38,  41,  42,  47,  48,  51,  55,  59].  The  dark 
solitons  appear  as  stable,  nondivergent  stripes  and  grids  in  the  transverse  plane  of  a  laser 
beam  that  passes  through  a  rectilinear  diffraction  screen  before  propagating  in  a  self- 
defocusing  nonlinear  material.  Materials  with  different  mechanisms  of  nonlinearity  manifest  the 
same  qualitative  results.  The  SDS  nature  of  the  observed  phenomenon  was  verified  by  numer¬ 
ical  simulations  of  the  (2+1  )-D  nonlinear  Schrodinger  equation,  analytical  solutions  for  the 
(1+1  )-D  case,  and  their  comparison  with  experimental  data.  This  work  was  done  by  this  Pi’s 
research  group  (including  G.  A.  Swartzlander  Jr.  and  H.  Yin)  in  collaboration  with  Dr.  D.  R. 
Andersen  of  the  University  of  Iowa  and  his  group. 

Dark  solitons  have  provoked  much  interest  since  they  were  first  shown  to  be  particular 
solutions  of  the  two-dimensional  (1+1  )-D  nonlinear  Schrodinger  equation  (NSE)  with  a  negative 
(self-defocusing  type)  nonlinear  coefficient  «2  (see  Eq.  (1)  below).  As  opposed  to  a  so-called 
"bright"  NSE  soliton,  which  propagates  as  a  stationary  wavepacket  of  finite  extent  in  a  (1+1  )-D 
nonlinear  medium  with  n2>0,  a  "dark"  soliton  is  characterized  as  a  stationary  "hole"  on  an  oth¬ 
erwise  uniform  plane  wave  --  it  exists  on  a  background  field  as  an  absence  of  energy  in  a  local¬ 
ized  region,  with  constant  size  and  shape  parameters,  and  remarkable  stability. 

So  far,  only  (1+1  )-D  temporal  dark  solitons  (i.  e.  intensity  minimums  propagating  along  a 
nonlinear  fiber  on  a  quasi-cw  bright  background)  have  been  observed  experimentally.  We 
observed  stable  spatial  structures  (e.  g.,  stripes,  crosses  and  grids)  in  the  transverse  cross- 
section  of  a  cw  optical  beam  propagating  in  a  material  with  a  self-defocusing  nonlinearity,  with 
these  structures  having  a  strongly  pronounced  soliton  nature  -  namely  that  of  spatial  dark  soli¬ 
tons  (SDS’s).  Although  no  (2+1  )-D  analytical  solution  for  dark  solitons  in  the  NSE  is  known  to 
date,  our  experimental  and  numerical  data  with  various  2-D  amplitude  and  phase  masks  pro¬ 
vided  strong  evidence  that  the  phenomenon  observed  by  us  is  indeed  due  to  spatial  dark  soli¬ 
tons.  Furthermore,  our  results  on  quasi-(1+1)-D  propagation  (see  also  [1])  have  shown  excel¬ 
lent  agreement  with  the  well  known  analytical  results  for  (1+1  )-D  dark  solitons. 

In  comparison  to  temporal  solitons,  SDS’s  are  easy  to  create  and  observe  experimentally, 
requiring  as  little  as  a  HeNe  laser  and  some  slightly  absorbing  fluid.  Various  applications  of 
SDS’s  can  be  envisioned,  such  as  optical  encoding,  limiting,  switching  and  computing,  and 


FIG.  1  Far-ficld  cross-sectional  images  for  propagation 
through  sodium  vapor  with  a  wire  mesh  in  the  incident  plane 
(column  1)  and  through  an  absorbing  liquid  with  the  “cross 
configuration  (column  2). 
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nonlinear  filtering. 

Our  exploration  into  SOS’s  formation  was  motivated  by  the  observation  of  intriguing  non¬ 
linear  transformations  of  far-field  (Fraunhofer)  diffraction  patterns  of  a  wire  mesh  placed  at  the 
input  face  of  a  sodium  vapor  cell  (see  Fig.  1(1  a)).  (Far-field  diffraction  is  highly  sensitive  to 
intensity-dependent  phase  changes  in  a  nonlinear  layer.)  In  all  the  cases  studied  so  far  with 
rectilinear  diffraction  screens,  the  linear  Fraunhofer  diffraction  pattern  evolves  into  various 
arrays  of  square  spots  as  the  laser  intensity  increases.  Subsequent  measurements  at  the  out¬ 
put  face  of  the  nonlinear  medium  (in  the  Fresnel  or  "near-field"  regime)  revealed  the  formation 
of  very  distinct  dark  stripes  which  had  caused  those  novel  nonlinear  far-field  patterns. 

Here  we  will  briefly  describe  our  earlier  results  and  then  reflect  in  more  details  more 
recent  results,  mostly  on  propagation  of  dark  solitons  inside  the  nonlinear  material  (near-field 
patterns),  which  also  included  not  only  "amplitude"  masks  (i.  e.  wire  mesh),  but  also  "phase" 
masks  (phase  step)  that  enabled  us  to  excite  most  fundamental  single  stripe  and  cross  dark 
solitons.  In  our  initial  experiment,  a  cw  frequency  stabilized  dye  laser  beam  was  passed 
through  a  wire  mesh  and  then  imaged  into  an  L=18/nm  long  cell  containing  sodium  vapor,  a 
strongly  nonlinear  material.  The  laser  frequency  was  tuned  slightly  below  the  D2  atomic  reso¬ 
nance,  corresponding  to  a  negative  nonlinearity  (n2  <  0).  The  Gaussian  beam  was  dissected 
first  into  a  -  3x3  array  of  spots  (see  Fig.  1  (/)).  The  spectacular  transformations  of  the  far-field 
patterns  described  above  were  observed  when  the  laser  was  scanned  from  an  off-resonance 
frequency  below  the  D2  transition  frequency,  toward  the  resonance  (thereby  varying  the  value 
of  n2)-  Fraunhofer  diffraction  was  observed  (see  Fig.  1(7, a)  far  off-resonance,  where  n2=0. 
However,  as  the  laser  frequency  approached  the  D2  line  (from  below),  the  far-field  profile 
transformed  itself  into  an  amazingly  well-organized  square  array  of  spots,  filling  the  central 
area  of  the  beam.  Fig.  1(7,b).  (At  the  self-focusing  side  of  the  D2  line,  where  n2>0,  the 
transformation  was  considerably  different  than  that  described  above  and  consistent  with  the 
earlier  results  for  self-focusing  materials.) 

The  geometric  beauty  of  the  far-field  patterns  and  their  stability  over  a  relatively  large 
range  of  intensities  and  driving  field  frequencies  led  us  to  believe  that  this  phenomenon  is  not 
attributed  to  the  specific  physics  of  the  nonlinearity  in  sodium  vapor,  but  rather  to  the  simple 
fact  that  the  nonlinear  component  of  refractive  index  is  negative,  i.  e.  n2  <0-  To  verify  this,  we 
tried  an  experiment  using  another  phenomenon  resulting  in  large  values  of  «2  <  0:  the  so- 
called  thermal  nonlinearity,  which  can  readily  be  induced  using  low  power  radiation  in  many 
slightly  absorptive  liquids.  The  results  of  these  experiments  showed  nonlinear  far-field  patterns 
amazingly  similar  to  each  other  and  to  those  of  the  sodium  vapor  experiment.  Although  the 
nonlinearity  due  to  the  thermal  effect  exhibits  some  spatial  nonlocality,  this  has  not  appreciably 
affected  the  observed  phenomenon,  presumably  because  the  characteristic  scale  of  the  nonlo¬ 
cality  was  smaller  than  the  soliton  size. 

The  simplest  and  probably  most  fundamental  wire  mesh  configuration  is  a  single  opaque 
"cross"  (see  Fig.  1(2)).  In  this  case  the  Fraunhofer  pattern  (Fig.  1(2,a))  experiences  a  nonlinear 
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transformation  into  a  grid  pattern  (Fig.  1  (2,b))  that  has  essentially  the  same  characteristics  as 
the  3  X  3  mesh  case  (except  that  there  are  fewer  spots),  indicating  similar  nonlinear  transfor¬ 
mations  in  both  cases.  To  understand  the  observed  phenomenon,  we  modeled  the  experiment 
using  the  (2+1  )-D  NSE  for  laser  beam  propagation  in  a  nonlinear  medium: 


2ik  dE/dz  +  VfE  +  k^n2 


E 


^  E/no  =  0 


(1) 


where  Vf  =d^/dx^  +d‘^/dy^  {  =  d^/dx^  for  the  (1+1  )-D  NSE),  and  y  are  the  transverse  coordi¬ 
nates,  z  is  the  propagation  coordinate,  E  is  the  complex  electric  field,  and  k  is  the  wavevector  in 
the  material.  Eq.  (1)  was  solved  by  us  using  numerical  methods.  The  simulated  far-field 
results  (see  Fig.  1(c))  show  the  same  features  as  the  experimental  far-field  patterns  in  Fig.  1 
(b).  The  agreement  between  the  experimental  and  numerical  results  is  remarkable. 

In  our  later  research  [1 ,4,24],  to  identify  the  physical  phenomenon  which  gave  rise  to  the 
observed  far-field  patterns,  we  studied  the  specific  features  of  wave  propagation  inside  the 
nonlinear  material.  This  was  accomplished  by  examining  the  near-field  patterns  at  the  output 
face  of  the  nonlinear  material  for  different  lengths  of  the  material,  and  for  various  boundary 
conditions,  including  a  single  wire,  two  parallel  wires,  two  orthogonal  sets  of  parallel  wires,  a 
wire  mesh,  a  single  phase  jump,  multiple  phase  jumps  intersecting  at  a  point,  and  two  parallel 
phase  jumps.  The  formation  of  pronounced  dark  stripes  or  grids  was  a  universal  phenomenon 
for  all  the  cases  studied.  In  general,  the  width  of  each  stripe  remained  almost  constant  as  the 
thickness  of  the  nonlinear  material  increased,  but  decreased  as  the  laser  field  strength 
increased.  The  number  of  these  dark  stripes,  which  tend  to  appear  in  pairs,  also  remained  con¬ 
stant  with  propagation  distance,  even  after  collisions.  These  observations,  together  with  the 
fact  that  n2<0  suggest  that  the  dark  stripes  are  spatial  dark  solitons.  Our  investigation  also 
shows  that,  notwithstanding  orthogonal  interactions,  (2+1  )-D  dark  soliton  stripes  behave  amaz¬ 
ingly  similar  to  the  analytical  (1+1  )-D  dark  solitons,  i.  e.,  it  appears  that  soliton  stripes  orthogo¬ 
nal  to  each  other  in  a  cross-sectional  plane,  propagate  almost  independently  of  each  other. 

To  verify  that  the  observed  phenomenon  was  indeed  attributed  to  (2+1  )-D  dark  solitons, 
we  investigated  the  two  most  fundamental  cases:  (/)  an  opaque  cross  ("amplitude  mask",  AM) 
composed  of  two  orthogonal  wires,  and  (li)  crossed  phase  steps  ("phase  mask",  PM)  con¬ 
structed  with  two  microscope  cover  slips  (see  Fig.  2  (1)  and  (2)).  The  corresponding  near-field 
images  are  shown  in  Fig.  2  for  a  thermal  nonlinear  material.  In  the  AM  case,  the  linear  Fresnel 
diffraction  pattern  (Fig.  2(1, a))  exhibits  a  gray  shadow  of  the  cross  flanked  by  bright  stripes.  In 
contrast,  the  shadow  is  completely  missing  in  the  experimental  nonlinear  profile  (Fig.  2(b)); 
instead,  two  high  contrast  dark  stripes  separated  by  a  distinct  bright  region,  are  formed  parallel 
to  both  axes  of  the  cross.  In  spite  of  the  magnification  in  size  with  increasing  beam  power  due 
to  self-defocusing,  the  width  of  the  dark  stripes  actually  decreases,  as  expected  for  dark  soli¬ 
tons. 


In  the  PM  case,  the  linear  Fresnel  diffraction  pattern  displays  a  broad  dark  cross  aligned 
with  the  axes  of  the  phase  steps,  flanked  by  diffractive  ringing.  As  the  laser  power  is  increased 
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(see  Fig.  2(2, b)),  the  width  of  the  central  cross  decreases  (even  at  the  intersection)  and  does 
not  split  into  two  (as  did  the  shadow  of  the  dark  opaque  cross  in  Fig.  2(b)).  This  also  occurs  in 
a  (1+1  )-D  nonlinear  system  when  the  boundary  condition  is  a  tc  phase  step,  with  the  central 
dark  stripe  referred  to  as  a  "fundamental  dark  soliton";  similarly,  the  pattern  in  Fig.  2(2b)  can  be 
regarded  as  a  ’’fundamental  dark  soliton  cross."  Our  numerical  solutions  for  both  nonlinear 
cases,  shown  in  Fig.  2(c),  again  reaffirm  that  Eq.  (1)  correctly  predicts  the  observed 
phenomena. 

A  strong  verification  of  the  soliton  nature  of  the  the  observed  phenomenon  was  achieved 
by  comparing  [1]  experimentally  measured  parameters  with  the  theory.  We  measured  the  soli¬ 
ton  divergence  angle,  6,  with  respect  to  the  optical  axis  (see  inset  of  Fig.  3),  to  determine  the 
soliton  characteristic  parameter,  A.„/  =  (the  soliton  amplitude  and  width  also  depend  on 

^„/),  where  ti„;  =  \n2  \  \E\^/2nQ.  When  a  single  wire  of  diameter,  is  small  compared  to 
the  beam  size,  then  's  determined  by  the  equation  X^i  =  cos(?i„/fc4)  ,  where  A  = 

The  excellent  agreement  between  the  data  and  theory  is  shown  in  Fig.  3.  The  inset  (Fig.  3) 
illustrates  the  formation  of  a  pair  of  diverging  dark  solitons  for  the  case  of  a  single  wire.  To  test 
the  integrity  of  dark  solitons  on  a  finite  (e.  g.,  Gaussian)  background,  we  computed  the  so- 
called  "soliton  constant" ,  and  found  that  this  parameter  was  nearly  constant  (=1 .76^//io«2^o)- 
This  indicated  that  these  dark  stripes  propagate  as  robust  formations,  even  though  the  back¬ 
ground  intensity  relaxes  adiabatically  due  to  self-defocusing  and  linear  diffraction. 

The  far-field  transformation  of  a  linear  Fraunhofer  pattern  into  a  tightly  organized  and 
ordered  nonlinear  pattern  (as  in  Fig.  1)  can  now  be  explained,  in  terms  of  the  formation  of  spa¬ 
tial  dark  soliton  stripes  and  grids.  This  effect  occurs,  in  general,  because  some  spatial- 
frequency  components  of  the  incident  beam  are  channeled  into  the  formation  of  solitons,  and 
thus  are  not  allowed  to  "radiate"  away  from  the  optical  axis  as  in  the  linear  case. 

In  conclusion,  we  have  observed  the  nonlinear  transformation  of  various  Fresnel  and 
Fraunhofer  diffraction  patterns  of  a  laser  beam  passing  through  a  rectilinear  amplitude  or 
phase  mask,  followed  by  a  self-defocusing  material.  Under  nonlinear  propagation,  spatial  dark 
solitons  formed,  appearing  as  dark  stripes  or  grids  in  the  beam  cross-section.  The  (2+1  )-D  sol¬ 
iton  structures  behaved  as  if  they  consisted  of  two  almost  independent  and  noninteracting 
(1+1  )-D  soliton  substructures.  The  similarity  between  our  experimental  data  and  numerical 
solutions  of  the  (2+1  )-D  NSE  was  remarkable. 

Amplification,  Instability,  and  Chaos  in  Nonlinear 
Counterpropagating  Waves 

In  the  pilot  paper  under  previous  AFOSR  support  and  mostly  in  the  extensive  exploration 
effort  under  this  grant  [3,5,22,32,40],  this  principal  investigator  and  Dr.  Chiu  Law  demonstrated 
for  the  first  time  that  two  linearly  polarized  counterpropagating  waves  in  a  Kerr  nonlinear 
medium  with  linear  dispersion  can  exhibit  amplification  and  multi-mode  temporal  instability. 


which  are  attributed  to  the  combined  effect  of  nonlinear  index  grating  and  linear  dispersion. 
We  also  explored  possible  use  of  this  effect  for  large  broad-band  amplification  that  may  exist 
in  the  system  (in  particular  in  a  nonlinear  optical  fiber)  if  the  pumping  is  below  the  threshold  of 
instability.  This  amplification  can  have  a  significant  potential  for  such  applications  as  optical 
fiber  communication. 

The  cross-interaction  of  two  counterpropagating  laser  beams  in  a  third  order  nonlinear 
material  is  a  conceptionally  simple  and  fundamental  process  in  nonlinear  optics.  In  the 
research  done  under  this  grant,  it  has  been  shown  [3,5,22]  that  regular  linear  dispersion  (i.  e. 
frequency  dependence  of  refractive  index)  can  be  a  natural  and  universal  agent  for 
amplification  of  perturbations  and  temporal  instability  in  the  counterpropagating  nonlinear 
waves.  It  is  well  known  that  linear  dispersion  can  give  rise  to  nonlinear  optical  effects  such  as 
spatial  instability  of  single  plane  wave,  formation  of  soliton  in  nonlinear  fiber,  and  amplification 
of  counterpropagating  waves.  This  suggests  that  linear  dispersion  in  combination  with  non¬ 
linear  processes  can  dramatically  change  the  dynamical  behavior  of  a  system. 

In  our  pilot  research  we  started  investigation  of  stability  of  eigenpolarization.  The  total 
complex  electric  field  of  the  two  counterpropagating  plane  waves  in  the  Kerr  nonlinear  medium 
was  represented  as:  [£1  (z,r)e'*^ -i- £2(2,0^"'^]  where  E^{z,t)  and  E2{z,t)  are 

the  slowly  varying  envelopes  of  forward  (+z)  and  backward  (-z)  propagating  waves  respec¬ 
tively.  The  dynamics  of  these  envelopes  is  governed  by  two  coupled  nonlinear  Schrodinger 
equations: 


•j.1 


-h 


dEj 


dt 


7  =  1.2 


(1) 


where  Ij  =  \Ej\^  is  the  intensity  of  the  respective  propagating  wave,  |i  =  is  the  linear 

dispersion  parameter,  P  =  n2k/n  is  the  nonlinear  parameter,  n  is  the  linear  refractive  index  at 
the  frequency  (o,  is  the  linear  group  velocity,  ^2  'S  the  nonlinear  refractive  index  coefficient, 
and  k  is  the  wave  number  in  the  medium.  The  coefficient  2  in  the  right-hand  side  in  Eq.  (1) 
reflects  light-induced  nonreciprocity. 

To  analyze  the  small  perturbation  stability  of  the  steady  state  and  also  the  amplification, 
we  represented  both  waves  in  the  system  as  slightly  perturbed  steady  state  in  the  form: 

Ej  =  £;o©[1  +  A  (2) 

where  y=1 ,2,  Ejq  is  the  steady  solution  for  Eq.  (1),  A  and  A2,;  are  normalized  amplitude  of 
the  small  perturbations,  'i  =  tVg/L  is  the  normalized  time,  and  ^=z/L  is  the  normalized  distance 
of  propagation,  with  L  being  the  total  length  of  nonlinear  material.  Substituting  Eq.  (2)  into  Eq. 
(1)  and  linearizing  the  latter  equation,  we  obtained  [3,5,22]  linear  propagation  equations  ior  Aj^j 
and  solved  them.  We  assumed  for  simplicity  equal  intensities  of  the  waves,  i.  e.  /io=^20=^- 
Using  the  boundary  conditions  for  the  perturbation  amplitudes  we  obtained  [3,5,22]  an  equation 


for  X.  and  solved  it. 

The  details  of  our  results  can  be  found  in  [3,5,22];  some  preliminary  new  results  have 
been  reported  by  us  previously.  We  will  briefly  summarize  our  results  on  instabilities  and  will 
give  a  bit  more  details  on  the  most  recent  results  on  amplification  (especially  in  fiber 
waveguides  pumped  by  the  counterpropagating  waves).  We  found  the  boundaries  of  instability 
in  the  space  of  parameters  of  the  system  (dispersion,  nonlinearity,  pumping  intensity  and  the 
length  of  the  nonlinear  system).  For  further  discussion  we  introduce  dimensionless  parameters: 
p  =  p/L  --  normalized  pumping  intensity  (note  that  p  can  have  negative  sign  depending  on  non¬ 
linearity  p),  and  d  =  [ivp{2L)  -  normalized  dispersion.  We  found  that  the  entire  domain  of  ins¬ 
tability  consists  of  multitude  of  individual  unstable  modes,  and  we  found  solutions  for  them. 
These  individual  modes  essentially  are  longitudinal  modes  in  a  light-induced,  distributed- 
feedback  resonator.  We  found  also  the  boundary  encompassing  all  the  unstable  solutions. 
The  (encompassing)  threshold  intensity  increases  as  dispersion  decreases;  our  numerical  cal¬ 
culations  show  that  for  sufficiently  small  \d\,  this  dependence  can  be  best  described  by  a 
surprisingly  simple  relationship 

Per  =  -sgn  ipd)  1  I  ri  I  I  (3) 

where  r)  is  numerically  determined  to  be  10.0+1%,  and  sgn(pri)  is  the  sign  of  pd.  The  neces¬ 
sary  condition  for  initiating  instability  (i.  e.  to  have  per  >  0)  is  that  the  signs  of  nonlinearity  nz 
and  dispersion  |x  must  be  opposite,  which  coincides  with  the  necessary  condition  for  formation 
of  a  soliton  and  spatial  instability  in  single-wave  propagation. 

When  the  pumping  is  only  slightly  above  the  threshold,  only  one  mode  become  unstable, 
then  two,  and  so  on.  When  the  number  of  modes  is  greater  than  one,  they  compete  with  each 
other  and  some  of  them  can  become  dominant  while  others  fade  out.  If  only  one  mode  is 
unstable,  it  eventually  develops  into  periodic  self-sustained  oscillations  with  a  stable  amplitude. 
When  the  value  of  \p/pcr  I  's  increased  above  1.28,  the  instability  and  resulting  self-sustained 
oscillations  develop  much  faster;  as  the  oscillations  develop  in  time,  the  second  subharmonic  is 
excited  which  later  evolves  into  aperiodic  oscillations  with  randomly  modulated  amplitude  and 
phase  indicating  onset  of  chaos.  This  behavior  further  develops  into  strongly  pronounced 
chaos  when  the  pumping  intensity  is  significantly  exceeds  the  threshold  of  instability,  Eq.  (3). 

To  estimate  the  threshold  intensity  for  instability,  consider  an  example  of  a  1  km  long 
single-mode  fiber  with  Ge-doped  silica  core  at  wavelength  1 .55  pm  with  group  velocity  disper¬ 
sion  D(k)  =  6.5x10“^  [note  that  \i=-D  {k)/{kc^)]  (the  corresponding  d  = -8.02x10"''^),  refrac¬ 
tive  index  n  =  1.44,  and  ^2  =3-2x10"^®  cm^/W  ,  in  the  lossless  approximation.  We  find  the 
threshold  intensity  4^,  Eq.  (4),  for  such  a  fiber  to  be  10^  W/cm^,  which  is  below  the  damage 
threshold  10^°  W/cm^  for  fused  silica.  Crude  estimate  shows  that  the  losses  existing  in  the 
real  fiber  (-0.5  dB/km  )  would  require  only  about  two  times  higher  threshold  intensity.  If  we 
use  the  SF-59  glass  with  ^2  =7x10“^®  cm^AV  a\  wavelength  1.06  pm,  and  assume  that  112  at 
wavelength  1 .55  pm  is  of  the  same  order  as  that  of  1 .06  pm  and  that  dispersion  is  roughly  the 
same  as  for  plain  glass,  the  critical  intensity  is  reduced  to  ~A25 kW/cm^  for  the  same  length. 


Figure  4 .  Pump-probe  configuration  for 
amplification  measurement  where  E\  and  E2  ure 
pumping  waves  and  input  probe 

wave  with  ^10^2.1,  E2oAi^2  and  ^20^2.2  generated 
as  a  result  of  the  wave  mixing  process. 


sn/n 


Figure  5.  Amplification  Gi  of  the  probe  wave 
versus  frequency  modulation  index  5Q/CI  for  the 
case  of  a  1  km  long  Ge-doped  silica  fiber.  (Curves: 
1  —  60%  of  instability  threshold  p„  [120  times  of 
the  amplification  threshold,  /««?]»  2  —  40%  of  pg^ 
[SOx  of  fainp])  3  ■—  20%  of  Pgf  [40  X  of  •/amp]}  und  4  “ 
5%  of  Pgr  [lO  X  of  /«mp]). 
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The  most  interesting  opportunity  provided  by  systems  in  consideration,  is  amplification  of 
a  weak  probe  signal  under  counterpropagating  pumping,  when  the  pumping  is  way  below  the 
instability  threshold.  To  find  the  gain  spectrum  of  this  system,  we  inject  a  weak  probe  beam 
into  the  nonlinear  medium  and  scan  the  beam  frequency.  In  Fig.  4,  we  show  a  configuration 
with  two  pump  beams  and  four  weak  beams.  Suppose  that  the  probe  beam  (£ioA  ^ )  has  fre¬ 
quency  of  Q-I-5Q,  i.  e.  its  frequency  deviated  from  that  of  the  pump  beams  with  normalized 
frequency  Q  by  -I-6Q:  then  we  expect  a  phase  conjugate  signal  (£^20^2,2)  fo  be  reflected  back 
with  frequency  Q-6Q. 

Our  numerical  results  are  shown  in  Fig.  5.  We  define  the  intensity  amplification  Gi  for 
the  probe  wave  A 1-1,  as:  Gi  =  |Aij(^=1)/Aij(^=0)|^.  The  variation  of  G 1  versus  normal¬ 
ized  frequency  detuning  50  is  depicted  in  Figs.  5  for  the  1  km  long  Ge-doped  silica  fiber  with 
pumping  at  60%,  40%,  20%,  and  5%  of  the  threshold  intensity  Icr  =  9.3MW/cm^.  One  can  see 
that  even  when  the  pumping  is  well  below  the  threshold  of  instability,  the  bandwidth  and 
amplification  are  still  appreciable.  For  instance,  in  the  case  of  pumping  at  20%  of  the  threshold 
of  instability  the  gain  =20  and  bandwidth  =2.5x10®  (7.4x10"'''  Hz  or  0.4%  of  the  pumping  fre¬ 
quency). 

Under  realistic  conditions  d<t:^ ,  we  found  a  good  analytical  approximation  for  the  gain  Gi 
[3,5,22].  We  shoved  that  the  gain  spectrum  differs  significantly  between  detuning  domains 
separated  by  the  detuning  {6Q.)o  =  {8\pd\)~‘'^  ln(1+/?^).  The  portion  of  the  gain  spectrum  for 
6Q<t:(5i^)o  corresponds  to  the  contribution  from  dispersionless  four  wave  mixing.  Beyond 
(5n)o,  the  variation  of  gain  obeys  -exp (25Q[2|rf/7  |-(cf6Q)^]^^)  when  pumping  power  \p  \  is 
substantially  below  pcr-  This  part  of  the  gain  is  solely  caused  by  dispersion-related  process. 

We  found  that  a  nonzero  amplification  (or  G-i  >  0)  can  be  obtained  for  any  pumping  inten¬ 
sity,  however,  for  each  fixed  intensity  |/7  | ,  the  frequency  detuning  should  be  smaller  than  cer¬ 
tain  cutoff  frequency,  (5Q)cr, 

|5Qi  <(5a)^,  =  V2|/7/rf|  . 

Therefore,  a  rough  estimate  of  the  amplification  bandwidth  is  given  by  (5Q)cr,  Eq.  (4),  which  is 
valid  for  arbitrary  pumping,  |p  |  < Pcr-  The  maximum  bandwidth  (6Q)max  ^  dispersion 
d  \s  attained  when  \p  \  =Pcr,  Eq.  (3),  i.  e. 

mUax  =  ^2\log^\d\/d\  . 

The  location  of  the  peak  gain  {?>Q.)opt’  '•  s-  ^he  frequency  of  the  mode  with  lowest  threshold  of 
instability,  can  be  approximately  estimated  as: 

The  results  are  qualitatively  close  to  the  exact  values  from  numerical  calculation  shown  in  Fig. 
5.  All  of  them  indicate  that  the  nonlinear  fiber  with  dispersion,  pumped  by  counterpropagating 
waves  has  great  potential  as  an  all-optical  amplifier  operating  in  cw  or  quasi-cw  regime  which 
may  find  application  in  optical  gyroscopes,  optical  fiber  communications,  etc.  One  possible 


limiting  factor  in  amplification  is  the  losses  in  the  fiber  which  have  not  been  included  in  our  cal¬ 
culation.  To  estimate  the  effect  of  losses  in  the  optical  fiber,  we  introduce  a  threshold  intensity 
of  amplification  lamp,  which  satisfies  the  condition:  aL  =  T{Iamp),  where  a  is  the  loss  coefficient 
and  T{Iamp)  =  \'^[G^{Iamp)]  's  the  growth  rate  at  the  amplification  threshold  in  the  lossless 
approximation,  i.  e.  the  gain  must  be  equal  to  loss  at  lamp-  We  obtained  [3,5,22]  a  remarkably 
simple  expression  for4;„p  at  {hQ.)opt'- 

/amp=a/(2n2^)  (4) 

which  does  not  depend  on  either  length  or  dispersion.  For  the  silica  fiber  with 
a=0.5dB/km=QA2km~^  and  n2k  =  ^  km~\  Hence,  the  threshold  of  amplification  is 

46.2ik Wcm^  (or  0.5%  of  the  threshold  of  instability  in  the  example  with  L  =  Man  and 
\L=^4psec/nm-km).  For  a  typical  fiber  with  effective  area  of  50  \im^  ,  the  threshold  pumping 
power  for  amplification  is  23mW.  This  excellent  projected  performance  is  due  to  the  low  loss  in 
the  fiber. 

To  conclude  this  section,  we  discovered  that  linear  frequency  dispersion  together  with 
Kerr  nonlinearity  can  result  in  amplification  and  temporal  instability  of  nonlinear  counterpro- 
pagating  waves.  As  the  pumping  increases  above  certain  threshold,  self-oscillations  are 
excited;  upon  further  increase  of  pumping  they  evolve  into  subharmonics  and  chaos.  With 
under-threshold  pumping,  large  gain  and  broad-band  amplification  are  to  be  found.  The 
mechanism  of  the  entire  phenomenon’ can  be  explained  in  terms  of  positive  distributed  feed¬ 
back  from  the  nonlinear  index  grating  formed  by  the  two  laser  beams.  Our  calculations  show 
that  the  amplifiers  based  on  the  nonlinear  optical  fiber  pumped  by  counterpropagating  waves 
with  relatively  low  power  have  great  potential  for  various  applications. 

2.iv.  Field-Gradient-Induced  Second  Harmonic  Generation  in  Vacuum 

In  this  research  we  revised  our  previous  result  on  QED-based  SHG  in  vacuum.  Our  new 
results  [11,15,20]  show  that  the  photon-photon  scattering  in  vacuum  can  give  rise  to  the 
second-harmonic  generation  of  intense  laser  radiation  in  a  dc  magnetic  field  (even  in  the  "box" 
diagram  approximation  of  quantum  electrodynamics  which  was  a  source  of  some  controversy 
in  the  literature)  if  the  symmetry  of  interaction  is  broken  by  nonuniformity  of  optical  wave  +  dc 
field  system  in  time  or/and  space.  Specific  examples  considered  by  us  were:  optical  pulse 
plane  wave  and  a  Gaussian  laser  beam  propagating  in  either  uniform  or  nonuniform  dc  field. 

Photon-photon  scattering  (PPS)  is  perhaps  one  of  the  most  fundamental  quantum  electro¬ 
dynamics  (QED)  processes  which  may  also  result  in  nonlinear  optical  effects  in  vacuum  such 
as  the  birefringence  of  the  refractive  index  seen  by  a  probe  field  under  the  action  of  either  a  dc 
magnetic  (or  electric)  field,  or  intense  laser  pumping,  multiwave  mixing,  and  merging  of  two 
photons  into  one  (i.  e.  sum  frequency  generation)  under  the  action  of  a  dc  field.  If  observed, 
these  effects  may  provide  a  fundamental  optical  test  of  QED.  All  of  these  effects  are  based  on 
the  lowest  order,  so  called  "box"  diagram  approximation.  Yet  the  required  optical  fields  are  still 
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enormously  high  and  not  presently  available.  It  has  also  been  apparent  that  a  do  field  (either 
electric  or  magnetic)  may  greatly  assist  the  interaction.  However,  perhaps  the  most  interesting 
dc  field-assisted  processes:  the  merging  of  two  photons  into  one  (e.  g.  the  second  harmonic 
generation,  SHG)  and  photon  "splitting"  (in  essence,  a  parametric  process)  in  the  presence  of  a 
dc  field  -  have  caused  a  long  standing  controversy  (see  review  in  [15,16]),  with  the  results  on 
the  probability  of  the  photon  splitting  and  merging  differing  by  many  orders  of  magnitude.  Most 
recently,  we  have  proposed  (in  the  work  under  the  earlier  AFOSR  grant,  Y.  J.  Ding  and  A.  E. 
Kaplan,  Phys.  Rev.  Lett.  63:  2725  (1989))  the  SHG  of  laser  radiation  in  vacuum  in  the  pres¬ 
ence  of  a  dc  magnetic  field.  As  result  of  the  discussion  (Phys.  Rev.  Lett.  2744-2746 
(1990))  following  that  publication,  we  came  to  the  conclusion  [11,b]  that  the  result  for  box 
diagram  vanishes  for  a  cw  plane  wave  and  an  uniform  dc  field  which  was  essentially  the  main 
case  considered  earlier;  we  also  suggested  that  the  nonvanishing  effect  may  result  from  the 
nonuniformity  of  dc  field. 

Thus,  aside  from  the  issue  of  immediate  observability,  the  fundamental  question  arises 
whether  the  vanishing  contribution  of  the  box  diagram  is  due  to  some  fundamental  laws  of 
QED,  or  only  to  the  chosen  configuration:  plane  monochromatic  (cw)  optical  wave  +  uniform 
dc  field.  In  our  most  recent  research  [11,15],  we  showed  that  the  nonvanishing  contribution  of 
the  box  diagram  can  result  from  the  nonuniformity  (or  gradient)  of  any  component  of  the  entire 
field  system  (in  particular  -  optical  field)  in  space  or/and  time.  In  particular,  we  considered  (i) 
plane  wave  modulated  in  time  by  a  pulse  with  an  arbitrary  profile  and  finite  duration,  and  (ii)  a 
cw  Gaussian  beam  in  magnetic  field  with  an  arbitrary  spatial  distribution,  and  show  that  in 
these  configurations  nonvanishing  SHG  in  the  lowest  (i.  e.  box)  approximation  can  result  from 
the  laser  beam  propagating  either  in  a  uniform  or  nonuniform  dc  fields.  The  total  number  of 
SHG  photons  in  such  a  system  is  - 10^®  - 10^"^  orders  of  magnitude  higher  than  that  due  to 
the  next,  hexagonal  diagram  contribution  [1 1,15],  such  that  SHG  appears  to  be  truly  of  the  "box 
diagram"  nature. 

The  Heisenberg-Euler  Lagrangian  for  PPS  can  be  expanded  as  L=Lq+L/^+Lq  +  •  •  ■  , 
where  Lo  =  {E^-B^)/2  is  a  linear  term,  and  L4  =  {^)[{E^ -B^f +  7 -^f]  -  first  non¬ 
linear  term  that  corresponds  to  the  box  Feynman  diagram.  Here 
^  =  a/45KBcr  =  2.6x^0~^^Gauss~^  is  a  nonlinear  interaction  constant  with  a=e^/hc  =  ^A37 
being  the  fine  structure  constant  and  Bi.r  =  moc^/efi=4.4x‘\O^^Gauss  -  the  QED  critical  field. 
Lq  corresponds  to  the  hexagonal  Feynman  diagram,  etc.  Following  standard  procedures  one 
obtains  the  macroscopic  equations  in  the  form  of  classical  Maxwell’s  equations: 
V  •  5  =  0,  V  ■  D  =  0,  and  V x^-f- (1/c)  a^/9r  =  0  V x^  - (1  /c)  aD/ar  =  0,  with  the  constitutive 
relations  between  the  electric  displacement  1d  and  magnetic  field  and  electric  field  ^  and 
magnetic  induction  B  in  the  form  D  and  =B  +1^^^ ,  where 

D^^  =  dL^/Bi='^(2a't+7bt),  =-dL^/bt='^{2a't -7bt),  (1) 

with  a=E^-B^  and  b  =^’~B.  A  cw  plane  wave  does  not  exhibit  any  nonlinear  effects,  since 
due  to  its  properties,  E^  =B^,  B  =0,  the  nonlinearity,  Eq.  (1),  vanishes.  This  "degeneracy" 


of  the  third-order  nonlinearity  can  be  broken  by  the  field  nonuniformity  that  can  give  rise  to 
second -order  nonlinear  optical  effects  in  the  presence  of  a  strong  static  field. 

In  our  most  recent  work  we  considered  both  spatial  [11, a],  [15]  and  temporal  nonunifor¬ 
mity  [15].  The  perhaps  simplest  and  fundamental  example  of  temporally  nonuniform  wave  is  a 
pulse  plane  wave  propagating  in  a  uniform  dc  magnetic  field  normal  (say  along  the  axis  y)  to  the 
wave  propagation  axis  {fide  =BQef).  Suppose  that  the  amplitude  (and/or  phase)  of  fundamen¬ 
tal  plane  wave  is  arbitrarily  modulated  in  time  with  a  (complex)  envelope  mi(\}/),  where 
is  a  "retarded  time",  and  polarized  parallel  to  the  dc  magnetic  field  (this  polari¬ 
zation  provides  maximum  effect).  The  plane  wave  solution  for  the  SHG  amplitude,  M2>  'S  found 
then  as 

'^2(V.>’)=-7(^)^1  BQe^{d[u^{\fi)']/d\]f}-y,  (2) 

assuming  that  the  field  Bq  \s  "turned  on"  at  y  =  0.  Eq.  (2)  demonstrates  the  same  dependence 
on  the  distance  of  interaction  (m2  “3')  as  for  SHG  in  a  "classical"  nonlinear  medium  with  ideal 
phase-matching,  with  the  significant  difference  being  that  SHG  is  proportional  now  to  the 
time -derivative  of  the  driving  envelope.  The  use  of  very  broad  spectrum  radiation  can  further 
enhance  the  SHG  effect. 

We  also  considered  spatial  nonuniformity  [11,  15]  of  nonplanar  (in  particular,  Gaussian) 
wave;  we  assumed  a  cw  wave  (it  is  clear,  however,  that  a  combined  time/space  nonuniformity 
may  significantly  enhance  the  nonlinear  interaction).  We  found  that  predominant  contribution 
to  SHG  is  due  to  the  transverse  component  of  dc  magnetic  field.  Briefly,  our  results  can  be 
summarized  as  following.  In  the  case  when  the  fundamental  wave  is  a  linearly  polarized  Gaus¬ 
sian  beam  [TEMqq  mode),  we  found  that  a  "forced"  analytical  solution  for  vacuum  SHG  con¬ 
sists  of  two  modes;  one  of  them,  the  SHG  main  mode,  is  similar  to  the  fundamental  (laser) 
Gaussian  beam,  while  another  one  is  a  higher  order,  so  called  "doughnut"  Gaussian  mode, 
with  zero  amplitude  at  the  center  of  the  SHG  beam. 

If  the  dc  magnetic  field  is  symmetric  with  respect  to  the  position  of  the  beam  waist  (y=0) 
and  extends  at  the  distance  much  greater  than  diffractional  distance,  we  showed  that  the  main 
SHG  Gaussian  component  vanishes  in  the  far  field  area,  which  is  consistent  with  the  conven¬ 
tional  nonlinear  optics.  Thus,  this  SHG  component  could  be  regarded  as  an  analog  of  "classi¬ 
cal"  SHG  process  whereas  the  doughnut  component  is  a  feature  peculiar  only  for  gradient  non¬ 
linearity. 

We  considered  a  few  specific  cases,  in  particular  uniform  and  nonuniform  dc  magnetic 
fields.  In  the  case  of  uniform  dc  magnetic  fields,  we  showed  that  in  most  of  the  cases  of 
interest,  the  doughnut  component  contribution  exceeds  that  of  the  main  Gaussian  component 
by  order  of  magnitude  or  more.  For  the  arrangement  whereby  the  waist  of  the  laser  beam  is 
located  exactly  in  the  center  of  the  dc  magnetic  field,  the  main  component  contribution  at  the 
exit  point  vanishes  completely.  In  this  case,  the  number  of  SHG  photons  peaks  at  the  waist  of 
the  beam;  however,  at  the  exit  point  the  SHG  output  decreases  and  tends  to  some  constant. 
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which  is  due  to  the  only  surviving  doughnut  component.  For  the  fixed  length  of  interaction,  the 
number  of  SH  photons  is  maximal,  when  the  waist  of  the  beam  is  positioned  at  either  the  out¬ 
put  or  input  point.  We  also  considered  a  nonuniform  dc  magnetic  field  in  the  form  of  a  magnetic 
dipole,  originated  by  two  thin  parallel  "magnetic"  wires  positioned  in  the  plane  orthogonal  to  the 
laser  beam.  We  showed  that  in  such  a  case,  the  SHG  out  can  be  significantly  enhanced. 

An  interpretation  of  nonvanishing  SHG  in  the  nonuniform  fields  is  that  the  nonuniformity 
allows  for  momentum  transfer  between  photons  and  dc  field  (which  would  ultimately  result  in 
the  recoil  of  material  system  generating  the  dc  field),  thus  breaking  the  symmetry  that  causes 
vanishing  interaction  of  completely  uniform  field  system.  This  explanation  could  be  directly  cor¬ 
roborated  e.  g.  by  direct  QED  calculations  of  SHG  by  two  collinear  photons  +  elementary 
source  (particle)  of  dc  field,  similarly  to  quasi-elastic  scattering  of  single  photon  at  a  Coulomb 
potential.  Examples  of  such  sources  could  be  protons  (or  heavy  nuclei)  or  neutrons  with  two 
collinear  photons  "SHG-scattered"  at  the  particle  spin  and  Coulomb  dc  (electric)  field.  In 
macroscopic  terms,  the  SHG  is  originated  by  an  elementary  multipole  source,  in  a  limited  spa¬ 
tial  volume  <*cXi ;  we  found  that  in  the  lowest  approximation,  the  source  is  a  dipole  for  a  spin 
and  a  quadrupole  for  a  Coulomb  field. 

Thus,  we  demonstrated  the  feasibility  of  field-gradient-induced  SHG  by  the  intense  laser 
radiation  in  a  dc  magnetic  field  in  vacuum;  the  effect  does  not  vanish  in  the  QED  box  diagram 
approximation  only  if  the  participating  fields  are  temporarily/spatially  nonuniform. 

2.V.  Nonlinear  Optical  Research  on  Semiconductor  Films 
and  Quantum  Well  Structures. 

Part  of  the  effort  of  this  PI  is  devoted  to  experimental  research  on  nonlinear  optics  of  sem¬ 
iconductor,  in  particular  thin  ZnSe  films  and  superlattices  and  quantum  well  of  various  composi¬ 
tions.  To  the  extend,  this  is  a  continuation  of  pioneering  studies  on  nonlinear  properties  of 
ZnSe  attributed  to  resonant  excitons.  Another  reason  for  this  interest  was  to  use  self-bending 
effect  (first  predicted  by  this  PI  in  1969)  for  the  first  time  for  nonlinear  spectroscopy  of  nonlinear 
refractive  index,  in  particular  in  ZnSe  (see  Section  2.v.1  below).  This  research  was  a  colla¬ 
borative  effort  between  research  groups  of  this  PI  and  Dr.  J.  Khurgin  at  the  same  department. 
The  main  thrust  of  the  rest  of  research  on  nonlinear  optics  of  semiconductors  was  concentrated 
on  asymmetric  quantum  wells  which  is  mostly  in  the  area  of  research  interests  of  Dr.  J.  Khurgin 
and  his  group;  it  will  be  addressed  here  only  very  briefly  (see  Section  2.v.2  below). 
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2.V.I.  Spectral  Measurement  of  Nonlinear  Refractive  Index  in  ZnSe 
Using  Self-Bending  of  a  Pulsed  Laser  Beam 

In  this  research  [2,21],  this  Pi’s  research  group  (including  Y.  J.  Ding,  G.  Swartzlander,  and 
C.  T.  Law)  collaborated  with  Dr.  J.  B.  Khurgin’s  research  group  on  first  spectral  measurements 
of  nonlinear  refractive  in  ZnSe  using  self-bending  of  a  pulsed  laser  beam.  The  effect  of  self¬ 
bending  has  been  predicted  by  this  PI  in  1969  and  subsequently  observed  experimentally  in 
pulse  regime  in  NaCI  crystal,  and  recently  in  CS2,  the  first  observation  and  measurements  in 
cw  regime  (using  sodium  vapor  as  a  nonlinear  medium)  has  been  done  by  G.  Swartzlander 
(then  a  grad  student)  and  this  PI  under  AFOSR  support.  One  of  the  natural  applications  of  that 
effect  was  to  use  it  for  direct  measurements  of  nonlinear  refractive  index.  Our  measurements 
of  the  nonlinear  refractive  index  (n2)  spectrum  of  ZnSe  near  the  band  gap  (Xg^;,  ~450  nm)  at 
77  K  showed  that  the  maximum  nonlinearity,  n2  ~  1 .9x  10“®cm^/W,  measured  by  us,  was 
anomalously  large  to  be  explained  by  conventional  thermally-induced  bandgap  shrinkage 
(TIBS). 

ZnSe  is  one  of  the  semiconductor  materials  with  large  nonlinear  refractive  index  especially 
in  the  blue  domain  (^.^^^  -450  nm  at  77  K).  Although  it  is  believed  that  this  nonlinearity  is  basi¬ 
cally  due  to  TIBS,  the  nonlinear  mechanisms  are  still  not  completely  understood.  The  meas¬ 
urements  of  the  spectrum  of  nonlinear  refractive  index  in  the  vicinity  of  resonant  transitions  (e. 
g.  excitonic  or  band  edge  transitions)  are  crucial  for  the  physics  and  device  applications  of 
these  materials.  Although  many  new  techniques  for  these  measurements  have  been  pro¬ 
posed,  none  of  them  was  simple  enough  to  be  used  for  measuring  the  entire  spectra  of  non¬ 
linear  refractive  index  in  semiconductors. 

The  self-bending  method  for  measuring  spectrum  n2(^)  has  several  advantages  over 
other  methods.  It  requires  in  principle  a  single  shot  per  data  point  (k)  (although  due  to  lack  of 
the  detection  equipment  required  for  that,  our  measurements  here  were  based  on  multi-shot 
averaging),  whereas  e.  g.  the  Z-scan  technique  inherently  involves  many  shots  and  requires 
high  mechanical  and  laser  shot-to-shot  in  which  two  (pump  and  probe)  laser  beams  are 
required  and  the  results  for  «2  are  obtained  indirectly  by  using  Kramers-Kronig  transformation. 
Finally,  the  sign  of  n2  can  be  immediately  and  directly  determined  from  the  direction  of  self¬ 
bending  effect. 

An  ideal  case  of  self-bending  occurs  when  a  slab  beam  with  a  "triangular"  spatial  intensity 
profile  in  the  beam  cross-section  propagates  through  a  Kerr-like  nonlinear  medium  with  the 
refractive  index  n=no+n2^.  where  hq  is  the  linear  refractive  index  and  n2  the  coefficient  of 
nonlinearity.  Since  a  nonlinear  prism  is  induced  in  the  beam  path,  the  beam  (in  the  transparent 
nonlinear  material)  will  be  self-deflected  in  the  far-field  region  by  the  angle 

(^) 

where  a  0  is  the  beam  size,  L  is  the  thickness  of  the  nonlinear  medium  and  Iq  is  the  laser  peak 
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Fig.  6  Tlie  measured  spectra  of  self-bending  angles  dfji,  (dots  -  experimental  data, 
solid  curve  -  the  least  square  fitting)  and  of  nonlinear  refractive  index  n2  (cir¬ 
cles  -  recalculated  experimental  data,  broken  curve  -  the  least  square  fitting) 
vs.  the  laser  wavelength.  Inset:  the  measured  spectra  of  absorption  (a),  curve 
1,  and  of  self-bending  angle  due  to  TIBS  (^nl)t>  curve  2,  for 
Io~  1.5x10^  W/cm^.  The  shaded  area  corresponds  to  strong  absorption  of 
the  laser  beam. 


intensity  penetrating  into  the  nonlinear  medium.  For  our  experiment  parameters,  our  calcula¬ 
tions  show  that  typical  value  for  the  relative  error  introduced  with  the  use  of  right-triangular 
model  instead  of  a  semi-Gaussian  profile  used  in  experiment  was  less  than  10%.  In  our  exper¬ 
iment,  we  observed  the  self-bending  effect  in  bulk  ZnSe  over  the  range  446  nm  -  460  nm  with 

—ft  0 

maximum  ^2  as  large  as  1 .9  x  1 0' 

The  details  of  experimental  set-up  and  the  parameters  of  the  laser  pulse,  lenses,  etc.,  can 
be  found  in  Refs.  [2,21].  Here  we  will  briefly  discuss  the  results  of  the  experiment.  The  meas¬ 
urement  of  the  spectrum  of  the  self-bending  angles  for  different  laser  wavelengths  is  shown  in 
Fig.  6  (the  dots  and  solid  curve).  For  this  measurement,  the  laser  peak  intensity  was 
/o~  1.1  X  ^Ql^W/cm^  and  X  ranged  from  446  nm  to  460  nm.  The  self-deflection  angle  9^^^. 
had  a  clear  maximum  at  -448  nm,  where  ~  3.5  mrad.  The  spectrum  n2(^)  can  be  deter¬ 
mined  using  Eq.  (1)  and  measured  spectrum  of  However,  our  measurement  of  absorption 
shows  that  the  absorption  may  become  significant  near  the  bandgap  (e.  g.  -  30%  at 

X,-448  nm,  see  inset  in  Fig.  6  and  text  below).  Because  of  that,  712  should  be  corrected  using 

L 

Eq.  (1)  where  now  the  term  LIq  is  to  be  replaced  by  ^l{x)dx=lQ  ( 1  -e““^)/a,  where  a(^)  is 

0 

the  absorption  coefficient  at  wavelength  X.  The  absorption  spectrum  of  laser  beam  in  ZnSe 
was  measured  at  peak  laser  intensity  /q- 1.5x10®  W/cm^  and  is  plotted  in  inset  of  Fig.  6. 
Based  on  this  procedure  and  using  original  measurements  of  self-bending  angles  (solid  curve 
in  Fig.  6),  we  had  calculated  the  specfrum  of  712  which  is  plotted  in  Fig.  6  (broken  curve).  We 
were  unable  to  measure  n2  for  the  laser  wavelength  below  -  446  nm  because  of  strong  band- 
edge  absorption  (shaded  area  in  Fig.  6).  Our  experimental  results  showed  that  ZnSe  materials 
demonstrated  Kerr-like  nonlinearity  with  no  saturation  up  to  the  maximum  laser  peak  intensity 
-  3.1  xIO®  W/cm^,  measured  in  the  experiment,  which  was  verified  by  the  fact  that  the  self¬ 
bending  angles  were  linear  in  the  laser  intensity,  Iq,  see  Eq.  (1). 

The  spectrum  of  nonlinear  refractive  index  obtained  by  us  could  be  immediately  used  in 
attempt  to  either  identify  the  mechanism  of  the  observed  effect  or  at  least  rule  out  some  of  the 
known  mechanisms.  Although  there  are  a  few  nonlinear  mechanisms  described  in  the  litera¬ 
ture,  the  commonly  believed  dominant  effect  contributing  to  nonlinear  refractive  index  near  the 
bandgap  is  TIBS.  Assuming  the  worst  case  scenario  whereby  no  heat  transfer  takes  place 
during  a  single  laser  pulse,  the  maximum  increase  of  the  local  temperature  in  the  sample  due 
to  heating  with  the  laser  pulse  energy  penetrating  into  the  nonlinear  medium  Ji^ser 's 

=  ( A7’)o  C  ,  ( A7’)o  =  J laser  O./  %  Cp  Uq  (2) 

o 

where  x  is  the  distance  from  the  entrance  interface  of  the  laser  beam  and  Cp  =  0.34  J/K  •  cm 
is  the  heat  capacity  of  ZnSe  at  80  K.  From  the  absorption  spectrum  shown  in  the  inset  in  Fig. 
6,  one  can  determine  that  the  maximum  absorption  coefficient  measured  over  446  nm  -  460 
nm  is  ~74c777“^  at  X~AAQnm.  In  our  experiment,  with  the  typical  7/a^g;.~2.2x  10"^  7oM/e, 
ao~100p.77i,  and  L~5Q0\m,  the  absorption  can  increase  the  temperature  of  the  sample  at 
the  boundary  {x  =  0),  Eq.  (2),  by  ( Ar)o  =  0.15  K  at  A,=446  nm.  Consequently,  the  band  gap 
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decreases  by  AEg  =  {dEg/dT)  ■  AT,  where  dEg/dT =-8x10“^  eV/K,  \.  e.  AEg  —  1 .2x10”^  eV. 
Using  Moss  rule  we  determine  the  contribution  to  the  spectrum  of  nonlinear  refractive  index 
due  to  TIBS  as  An  =  {- dEg/dT)  hq  AT /4Eg.  Using  now  Eq.  (2),  it  can  be  translated  into  the 
self-bending  angle  due  to  TIBS  as 

{QNL)T  =  Jiaserno{-dEg/dT)i1-e-^^^^^)/4nEgCpal  (3) 

where  the  linear  dispersion  of  refractive  index  is  negligible  here.  Spectrum  of  {Qml)t  calculated 
from  Eq.  (3)  where  no  =  2.8,  is  plotted  in  the  inset  in  Fig.  6.  The  peak  0^,^  measured  by  us  in 
Fig.  6  occurs  at  X~  448  nm  and  0;^^  reaches  3.5  mrad.  At  the  same  wavelength,  the  self¬ 
bending  angle  calculated  based  on  TIBS  is  only  {Qnl)t  ~  ^  .3x10“^  rad,  i.  e.  at  least  two  ord¬ 
ers  of  magnitude  smaller.  In  order  to  further  verify  that  the  spectrum  of  the  self-bending  angle 
measured  by  us  was  not  due  to  TIBS,  it  is  instructive  to  also  compare  the  shapes  of  the  spec¬ 
tra  of  both  mechanisms.  This  comparison  shows  that  whereas  the  spectrum  for  TIBS  plateaus 
in  the  vicinity  of  fundamental  bandgap  (see  inset  in  Fig.  6),  our  measured  spectrum  has  a  dis¬ 
tinct  maximum.  TIBS  due  to  repetitive  laser  pulses  can  also  be  ruled  out  since  far-field  inten¬ 
sity  profile  was  independent  of  repetition  rate  (1  Hz  - 10  Hz).  Thus,  we  can  rule  out  TIBS  as  a 
mechanism  of  the  observed  nonlinearity.  Saturation  (band  filling)  and  free-carrier  effects  can 
also  be  ruled  out  since  they  both  result  in  negative  «2-  Further  investigation  is  apparently 
required  to  determine  the  real  nature  of  the  observed  nonlinearity. 

In  conclusion,  the  nonlinear  refractive  index  spectrum  of  ZnSe  was  directly  measured 
using  the  self-bending  of  a  pulsed  laser  beam.  It  was  demonstrated  that  this  effect  can  provide 
a  simple  and  reliable  method  of  direct  measurement  of  the  nonlinear  spectra  of  semiconductor 
materials.  Self-bending  in  semiconductors  has  potential  to  be  used  for  a  new  type  of  efficient 
nonlinear  optical  power  limiters,  in  particular  for  radiation  protection  of  optical  sensors  and 
resonatorless  optical  bistability. 

2.V.2.  Excitonic  nonlinearities  in  narrow  asymmetric  coupled  quantum  wells. 

This  Pi’s  research  group  has  also  been  involved  in  the  research  related  to  excitonic  non- 
linearities  in  narrow  asymmetric  coupled  quantum  wells;  the  research  was  motivated  by  expec¬ 
tations  that  narrow  asymmetric  coupled  quantum  wells  as  part  of  SEED  devices  may  under 
certain  circumstances  exhibit  a  blue  Stark  shift  (instead  of  the  regular  red  one)  which  might  not 
rely  on  sharp  excitonic  transition  usually  required  for  the  red  Stark  shift.  It  was  proposed  that 
blue-Stark  shift  can  be  achieved  by  a  proper  design  of  asymmetric  quantum  wells,  whereby  the 
first  two  heavy-hole  states  in  two  coupled  wells  are  approaching  the  resonance.  In  the  experi¬ 
ment  [6],  a  blue  shift  of  the  apparent  heavy-hole  (HH)  excitonic  transition  peak  in  the  photo¬ 
current  spectra  has  been  observed,  for  the  first  time,  in  narrow  GoAs /Al Q4Ga q  qAs  asym¬ 
metric  coupled  quantum  wells  (ACQW’s)  near  the  HH  pre-resonance.  With  an  external  reverse 
bias  of  only  -  2.35  V,  a  maximum  upward  shift  of  the  transition  energy  of  ~  6.1  meV  has  been 
measured  at  78  K.  This  blue  shift  is  anomalously  large,  which  cannot  be  explained  by  the  pre- 
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anticrossing  of  HH  energy  levels  in  two  coupled  quantum  well.  Sharp  change  of  the  linewidth 
of  the  HH  excitonic  transitions  has  also  been  observed  due  to  the  fluctuations  in  QW’s  and  bar¬ 
rier  widths  resulting  in  inhomogeneous  broadening  of  the  transitions.  CW  optical  bistability  was 
obsen/ed  with  the  external  feedback  due  to  this  inhomogeneous  broadening. 

Furthermore,  the  attempt  was  made  [8]  to  use  all  the  advantages  of  the  blue  shift  by 
attaining  an  intrinsic  feedback  mechanism  instead  of  a  conventionally  used  external  resistor. 
This  would  result  in  better  potential  for  integration  and  fabrication  of  devices.  In  the  experiment 
[8],  strong  excitonic  nonlinearity  in  the  photoconductive  response  of  a  PIN  photodiode  incor¬ 
porating  narrow  asymmetric  coupled  quantum  wells  (ACQW’s)  has  been  observed  at  78  K. 
When  the  PIN  photodiode  is  over-biased,  the  heavy  hole  (HH)  energy  levels  in  the  two  coupled 
quantum  wells  are  brought  into  resonance  by  increasing  the  laser  intensity.  Also,  both  the  light 
hole  (LH)  and  HH  excitonic  transitions  undergo  intensity-dependent  shifts.  Both  of  these 
effects  indicate  intrinsic  change  of  bias  due  to  redistribution  of  photogenerated  carriers,  and 
therefore,  the  existence  of  an  intrinsic  feedback  mechanism.  The  magnitude  of  the  blue  shift  of 
the  HH  excitonic  transitions  significantly  increases  when  the  laser  intensity  was  changed  from 
9.2  mW/cm^  to  ~  270  mW/cm^. 


2.vi.  Other  Research 

This  principal  investigator  continues  to  be  active  in  a  few  areas  of  nonlinear  optics  to 
which  he  made  pioneering  contribution  (in  some  of  them,  like  nonlinear  interfaces,  sixteen 
years  ago).  Most  recently,  he  reviewed  the  most  recent  research  and  outlined  new  research 
directions  on  nonlinear  optics  of  a  single  slightly-relativistic  electron  [16]  (with  Y.  I.  Ding),  on 
bistable  solitons  and  their  applications  to  optical  switching  [13]  in  collaboration  with  R.  H.  Enns, 
S.  S.  Rangnekar  of  Simon  Fraser  University  (Canada),  and  on  nonlinear  interfaces  in  colla¬ 
boration  with  P.  Smith  and  J.  Tomlinson  of  Bellcor  [17,18]. 
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